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Abstract 



The first observation of the Z boson decaying to four leptons in proton-proton 
collisions is presented. The analyzed data set corresponds to an integrated lu- 
minosity of 5.02 fb -1 at y/s = 7TeV collected by the CMS detector at the Large 
Hadron Collider. A pronounced resonance peak, with a statistical significance of 
9.7 cr, is observed in the distribution of the invariant mass of four leptons (elec- 
trons and /or muons) with mass and width consistent with expectations for Z bo- 
son decays. The branching fraction and cross section reported here are defined by 
phase space restrictions on the leptons, namely, 80 < < 100 GeV, where mu 
is the invariant mass of the four leptons, and mm > 4 GeV for all pairs of leptons, 
where mu is the two-lepton invariant mass. The measured branching fraction is 
B(Z -4 U) = (4.2+[j|(stat.) ± 0.2(syst.)) x 10~ 6 and agrees with the standard model 
prediction of 4.45 x 10~ 6 . The measured cross section times branching fraction is 
(7(pp -> Z)B(Z -4 U) = 112+2o(stat.)t5(syst.)+2(lumi.)fb, also consistent with the 
standard model prediction of 120 fb. The four-lepton mass peak arising from Z — > M 
decays provides a calibration channel for the Higgs boson search in the H — > ZZ — > M 
decay mode. 
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1 Introduction 

We present the first observation of Z — > A£ decays in pp collisions at a/s = 7 TeV, and measure- 
ments of the production cross section times branching fraction for Z decays into four leptons 
c"(pp — > Z) B(Z — > 4c£) and the branching fraction B(Z — > 4£). In this paper, £ denotes an elec- 
tron or muon. Previously all four LEP collaborations reported observations of final states with 
four fermions, e + e~ — > Af, which include e + e — > Z — > 4/ (lHIl- However, the observation of 
Z — > A£ decays in pp collisions is of special interest. The clean resonant peak in the four-lepton 
invariant mass distribution at m±i = mz can be used for direct calibration of the four-lepton 
mass scale, the four-lepton mass resolution, and the overall four-lepton reconstruction effi- 
ciency in phase space similar to the light Higgs boson four-lepton decays, H — > ZZ — > 4£. The 
pp — > Z — > A£ process and its implications for the H — > ZZ — > A£ search at the Large Hadron 
Collider (LHC) were first studied in Ref. 0- 

These results are based on the data collected in 2010 and 2011 by the Compact Muon Solenoid 
(CMS) detector |6| at the LHC. The data set used in the analysis corresponds to a total inte- 
grated luminosity of 5.02 ± 0.11 fb _1 . 

2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid 13 m in length and 6 m 
in diameter providing an axial magnetic field of 3.8 T. Within the field volume are the silicon 
pixel and strip tracker, the crystal electromagnetic calorimeter (ECAL), and the brass/ scintillator 
hadron calorimeter (HCAL). Charged particle trajectories are measured by the tracker, cov- 
ering < (p < 2n in azimuth and \tj\ < 2.5, where the pseudorapidity is defined as tj = 
— ln[tan(0/2)], with 6 being the polar angle of the trajectory of the particle with respect to the 
anticlockwise beam direction. Muons are measured in the pseudorapidity range | tj \ < 2.4, with 
detection planes made using three technologies: Drift Tubes, Cathode Strip Chambers, and 
Resistive Plate Chambers. Matching the muons to the tracks measured in the silicon tracker 
results in a transverse momentum (pj) resolution of about 2% in the muon pj range relevant 
for the analysis presented in this paper (from a few to a few tens of GeV). The ECAL consists 
of nearly 76000 lead tungstate crystals distributed in the barrel region (EB, \tj\ < 1.479) and 
two endcap regions (EE, 1.479 < \t]\ < 3), and has an ultimate energy resolution better than 
0.5% for unconverted photons with transverse energies above 100 GeV. The electron energy 
resolution is 3% or better for the range of energies relevant for the measurement reported in 
this paper. A two-level trigger system selects the most interesting events for use in the offline 
physics analyses. A more detailed description of the CMS detector can be found in Ref. 0. 

3 Signal definition 

The leading-order (LO) Feynman diagram for the Z production and decay into four leptons, 
qq — » Z — > A£, is presented in Fig. [l] (left). This process is sometimes referred to as single- 
resonant four-lepton production. 

We define signal events as those with four leptons, e + e~e + e~ (4e), p. + p.~}i + }i~ (4f/), or e + e~}i + }i~ 
(2e2}i), with four-lepton invariant mass satisfying 80 < < 100 GeV and dilepton masses 
mil > 4 GeV for all six possible pairings of leptons. The lower limit on ma reduces background 
that rises rapidly as decreases. The branching fraction and cross section reported here are 
defined by these phase space restrictions. 
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4 Signal extraction 



Predicted partial widths and branching fractions for Z boson decays to 4e, 4]i, and 2e2}i final 
states are summarized in Table [Tj The results are obtained at LO using CalcHEP 3.2 Q, which 
takes quantum mechanical interferences into account. The partial width for the 2e2}i channel 
is different from twice the width in either the 4e or Ap. channel because decays to four leptons 
of the same flavour involve additional Feynman diagrams with permutations of same-sign 
leptons. We do not have a sufficient number of events to measure the differences between 
decay rates for the three four-lepton final states. We therefore measure the overall branching 
fraction B(Z — > A£) and assume the 4e, Ap, and 2e2p relative branching fractions predicted by 
theory (Table 1). The main irreducible background is the process qq — > Zy* — > A£, for which 
the LO Feynman diagram is shown in Fig. [TJ (right). Events are simulated with the next-to- 
leading-order generator POWHEG |8] . The effects of multiple pp collisions within each bunch 
crossing are taken into account in all simulated samples. 

Table 1: Partial widths and branching fractions for Z boson decays to 4e, Ap, 2e2p final states 
with nifj > A GeV for all lepton pairs. The branching fractions are calculated with CalcHEP 
3.2 at LO using the total Z boson width r to t = 2.4952 GeV [9|. Theoretical uncertainties are 
smaller than experimental uncertainties and are not shown in the table. 



Quantity of interest 


4e 


Ap 


2e2f/ 


U 


Partial width, T, (keV) 


2.95 


2.95 


5.21 


11.12 


Branching fractions, T; / Ttot 


1.18 x 1CT 6 


1.18 x 1CT 6 


2.09 x 10~ 6 


4.45 x 10" 6 


Relative fractions, f\ = YJY^g 


0.2655 


0.2655 


0.4690 





4 Signal extraction 

The trigger and selection criteria used in this analysis closely follow the H — > ZZ — > A£ search 
by CMS 1 10 1 . We use data collected with dielectron and dimuon triggers selecting events with 
at least two electrons or two muons with transverse momentum pj > 17 and 8 GeV. To match 
the trigger selection, we require that at least two leptons reconstructed offline have pj > 20 and 
10 GeV. In this phase-space region, we expect and observe a high trigger efficiency of 96-99%, 
depending on the final state. 

Muon candidates are reconstructed using two algorithms, one in which tracks in the silicon 
strip tracker are matched to hits in the muon detectors, and another in which a combined fit 
is performed to signals in both the silicon strip tracker and the muon system |TTf . The muon 

q i 1 t 




Figure 1: (Left) Diagram of the Z — > 41 process. (Right) Diagram of the Z7* —> 4£ process for 
the irreducible background of Z — > 2£ production with the initial-state radiation undergoing 
an internal conversion 7* — > 2£. Both Z and 7* are present in all propagators. The choice of 
propagators shown in the figures corresponds to the dominant contributions in the phase space 
80 < m u < 100 GeV. 
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candidates are required to be successfully reconstructed by both algorithms. Other identifica- 
tion criteria (the number of measurements in the tracker and muon systems, the fit quality of 
the muon track, and the consistency of the track with the primary vertex) are also imposed on 
the muon candidates to reduce the misidentification rate. The vertex with the largest value of 
Py for the associated tracks is chosen to be the primary vertex. According to simulation, this 
requirement provides the correct assignment for the primary vertex in more than 99% of both 
signal and background events. The average number of reconstructed vertices per event in 2011 
data is approximately 7. 

Electron reconstruction also involves two algorithms, one in which energy clusters in the elec- 
tromagnetic calorimeter are matched to signals in the silicon pixel detector and another in 
which tracks in the silicon strip tracker are matched to ECAL clusters |fl2"fl. The electron can- 
didates used in the analysis are reconstructed by either algorithm. Additional identification 
criteria (the spatial distribution of the shower in the ECAL, matching of the electron track with 
a cluster in the ECAL, and consistency of the track with the primary vertex) are imposed on the 
electron candidates to reduce the misidentification rate. Requirements are imposed on the ge- 
ometrical properties of the electron track with respect to the primary vertex and neighbouring 
tracks, as well as requirements on the number of missing hits in the innermost tracker layers in 
order to remove electrons produced by photon conversions in the detector material. 

Electrons are required to satisfy pj > 7 GeV and \tj\ < 2.5, while muons are required to satisfy 
p T > 5 GeV and \t]\ < 2.4. 

Leptons produced in the decay of Z (or 7*) bosons are typically isolated from hadronic activity 
in the event. The scalar sum of transverse energy deposits in the calorimeters and the transverse 
momenta of tracks in a cone of radius AR = a/ (A</>) 2 + (A?/) 2 = 0.3 in rj-<p space around each 
lepton is computed. This sum is corrected for the contribution from the other lepton candidates 
as well as for the average hadronic activity in an event. The ratio of this corrected sum to the 
lepton pj is required to be smaller than 0.275. 

To suppress background events with leptons originating from secondary vertices (such as 
heavy quark decays or misidentified pions), the significance of the impact parameter of each 
lepton relative to the primary vertex, S3D = £?3d/ C3D, is required to satisfy |s3d| < 4, where c^d 
is the impact parameter in three dimensions and cr 3D is its associated uncertainty. 

Events are selected with at least four leptons satisfying the above criteria. The same-flavour, 
opposite-charge lepton pair with mass closest to the nominal Z boson mass is identified. A 
second pair is constructed from the remaining leptons with the highest pj satisfying the re- 
quirement of the same flavour and opposite charge. We require all six dilepton combinations 
constructed from the four selected leptons to satisfy w% > 4 GeV and the four-lepton invariant 
mass must be in the range 80 < m 4 £ < 100 GeV. The overall reconstruction and selection effi- 
ciencies for simulated signal events are found to be 25%, 59%, and 33% for the 4e, 4/y and 2e2}i 
channels, respectively. The overall theoretical acceptance efficiencies (per channel, in the same 
order) are 8.9%, 14%, and 11.1%. 

We estimate the irreducible four-lepton background from simulation using selection efficiencies 
corrected for differences between the simulated efficiencies and those measured in data. 

Processes with one or more misidentified or nonprompt leptons contributing to the signal re- 
gion 80 < < 100 GeV are dominated by Z + X production. Their rate is measured in data 
using control regions and is found to be negligible. The methodology used is identical to that 
used in the H -> ZZ -> U search |Tl3|. 
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5 Results 



We use pp —> Z — > 2]i events as a control sample in order to cancel systematic uncertainties in 
the measurement of the branching fraction B[Z — > A£). The muon selection criteria are identical 
to those in the Z — > A£ event selection with pj thresholds pj > 20 and 10 GeV applied to the 
higher and lower pj muons, respectively. 

5 Results 

5.1 Observation of Z — > 4€ decays 

Figure [2] shows the four-lepton mass distribution for events that pass all the selection criteria 
except the invariant mass requirement. Events with all three final states, 4e, 4^/, and 2e2ji, are 
included. The prominent peak at = niz constitutes the first observation of pp — > Z — > A£. 
The number of events with 80 < < 100 GeV is 28, in good agreement with the standard 
model expectation of 30.0 events. Table|2]shows good agreement between the expected and ob- 
served events in each of the 4e, 4fi, and 2e2^ channels. The probability for a background fluc- 
tuation to be at least as large as the observed maximum excess can be evaluated by generating 
sets of simulated data incorporating all statistical and systematic uncertainties. This probabil- 
ity can be expressed as a corresponding number of standard deviations using the one-sided 
Gaussian tail convention. The results presented in this paper are obtained using asymptotic 
formulae [14] and include the systematic uncertainties described in the next subsection. The 
statistical significance of the peak in Fig. [2] is 9.7c. 

Table 2: Numbers of expected and observed events. The yields are determined by means of 
fits described in the text; the fit results in the 4e channel are not meaningful because only two 
events are selected. 



Final state channels 


4e 


4 F 


262^ 


U 


Irreducible background (pp — > Zj* — > A£) 


0.07 


0.25 


0.14 


0.46 ± 0.05 


Other (reducible) backgrounds 


0.01 


0.01 


0.05 


0.07 ±0.1 


Expected signal (pp — >■ Z — > A£) 


3.8 


13.6 


12.0 


29.4 ± 2.6 


Total expected (simulation) 


3.9 


13.9 


12.2 


30.0 ± 2.6 


Observed events 


2 


14 


12 


28 


Yield from fit to the observed mass distribution 




13.6 ±3.8 


11.5 ±3.1 


27.3 ± 5.4 



5.2 Measurement of <r(pp Z) B(Z U) and B(Z U) 

The numbers of Z — > A£ (signal) and Z — ?► 2]d (control) events passing the selection may be 
expressed as follows: 

N ie -N*l g = J £L<r(pp^Z)B(Z^4£)f i ef cc e f : xp c t (/ = 4e,4^,2e2 F ), (1) 

i 

N^-Nj* = La(pp^Z)B(Z^2 ¥ )ef^ 2}l e e ^ 2p c z ^ (2) 

where 

• L denotes the integrated luminosity; 

• o"(pp — > Z) = 26.9 nb is the theoretical Z boson production cross section (80 < < 
100 GeV), calculated with FEWZ El at NNLO level; 

• B{Z — > 4£) is the signal decay branching fraction (for m^> 4 GeV); 



5.2 Measurement of u(pp ->■ Z) B(Z ->■ 4£) and B(Z -> 4£) 
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• H(Z -> 2p) = 0.03366 ± 0.00002 is the Z -4 2p branching fraction 0; 

• /, is the relative fraction of all A£ events in the zth channel (i = 4e,4^, 2e2^) from 
Table [1} 

• ef cc is the theoretical acceptance of the lepton pj and r\ requirements used in the 
analysis; 

• e^ xp is the experimental efficiency, as obtained from the simulation, to reconstruct 
events within the acceptance; 

• c, is the data-to-simulation correction factor for the experimental efficiency derived 
from Monte Carlo events. The C{ are within one percent of unity. 

Equation |T]) allows the extraction of the production cross section times branching fraction 
cr(pp -4 Z) B(Z -4 U) = 112+^ (stat.) fb, while the ratio of Eqns. jl} and Q allows extrac- 
tion of the branching fraction H(Z —¥ Al) = (4.2lg| (stat.)) x 10 -6 with cancellation of several 
systematic uncertainties. 

The signal and the irreducible background have the following systematic uncertainties: 

• Interference between the signal and background production is less than 0.2% in the 
range 80 < < 100 GeV, and is accounted for as a source of systematic uncer- 




60 80 100 120 140 160 

m 4£ [GeV] 



Figure 2: Four-lepton invariant mass distribution for events passing all selection requirements 
except that on m^g. The data are shown by points. The filled histograms represent standard 
model expectations for pp — > Z/Z7* — > A£ and for reducible backgrounds. The three final 
states, 4e, 4]i, and 2e2ji, are combined. 
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5 Results 



tainty. 



• Uncertainties in parton distribution functions lead to uncertainties in the total cross 
section and acceptance. For pp — > Z7* — > 4£ and pp — > Z — > 2\i, these uncertainties 
are 4% [10] and 1% fll|. For pp ->■ Z -» 4£ we obtain an uncertainty of 4% by 
following the same procedure as in Ref. [fTUI . 

• The theoretical uncertainties related to the choice of QCD renormalization/ factorization 
scales are evaluated by varying these scales up and down by a factor of two. For 
pp — > Z7* — > 4£ and pp — > Z — > 2]i, these uncertainties are both 2% IfTUl [T6| . For 
pp — > Z — > 4£ we obtain an uncertainty of 1% by following the same procedure as 

in Ref. ffToH . 

• The effect of bremsstrahlung radiation on event acceptance is estimated to be negli- 
gible. 

• Lepton reconstruction, isolation, and impact parameter uncertainties as well as un- 
certainties associated with the high-level trigger efficiency are determined using a 
"tag-and-probe" method [16] where efficiencies and uncertainties are measured in 
samples of Z — > 2.1 events. The statistical uncertainties on the lepton efficiencies 
obtained in different (pj, f]) bins are propagated to four-lepton reconstruction un- 
certainties for Z — > 4£ and Z7* — > 4£ events. These uncertainties range from 1% to 
6%. Correlation between the "tag-and-probe" efficiency measurement at the Z — > 2i 
calibration signal are taken into account. 

• The uncertainty on the integrated luminosity is 2.2% fF7\. 

The reducible background is estimated directly from data. The predicted yield of the reducible 
backgrounds is approximately 0.1 event (less than one percent of the signal yield, as shown in 
Table [2) and has a total uncertainty, also evaluated from data, of 130%. 

To account for systematic uncertainties, including correlations among different channels, as 
well as between signal and background, we construct the likelihood for the four observed final 
states (4e, 4ji, 2e2ji, and 2]i). Using the full likelihood in conjunction with the profile likeli- 
hood method [14J, we measure the production cross section times branching fraction and the 
branching fraction for Z — > 4£ to be 



These measurements agree with the standard model predictions, 120 ± 5 fb and 4.45 x 10 6 , 
respectively. 

5.3 Z — )• 4€ decays as a calibration channel in the H — > ZZ — >■ M search 

The Z — > 4£ decays give a narrow resonant peak in the four-lepton invariant mass distribution, 
which can be used as a calibration channel in the context of the Higgs boson search in the 
H — > ZZ — > A£ decay mode. Figure [3] (left) shows that the number of events in the Z — > A£ peak 
at m^£ = mz is approximately an order of magnitude larger than the expected number of events 
for the SM Higgs boson with a mass wjh = 125 GeV [fl"3llT8| . The Z — > 4£ peak can be used 
for a direct calibration of the four-lepton mass scale, the four-lepton mass resolution, and the 
overall four-lepton reconstruction efficiency in a phase space similar to the Higgs boson four- 
lepton decays. Such a direct calibration using the Z — > 4£ peak is complementary to the current 



cr(pp^Z)6(Z^4£) 
B(Z -> U) 
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indirect tag-and-probe method making use of Z — > 11 events, with the benefit that Z — > A£ has 
similar kinematics to H — > ZZ — > 4£. 

Figure [3] (right) shows the results of an unbinned maximum-likelihood fit to the four-lepton 
mass distribution for the observed events. The fitted signal yield is shown in the last line of 
Table |2j along with the results of separate fits to the mass distributions in the three channels. 
The background shape comes from a simulation of pp — > ZZ — > A£ events with the overall 
normalization floating in the fit. Differences in shape between the reducible and irreducible 
backgrounds are ignored because of the small contribution of the former. The signal shape is 
a convolution of the Breit-Wigner and Crystal Ball lTl9l functions. The central value and width 
of the Breit-Wigner function are fixed at the Z boson mass and width [9]. The Crystal Ball 
parameters are free to vary in the fit. Data used for this analysis is only a (lower energy) subset 
of that used in the recent discovery paper |fl3| . 

From the fit results, as shown in Fig. [3] (right), one can see that the offset of the peak is 0.4 ± 
0.4 GeV (the "mean" value of the Crystal Ball function represents an offset with respect to the 
fixed Breit-Wigner function peak position), or, in relative units, 0.4 ± 0.4%. These numbers can 
be used to constrain the possible systematic uncertainty of the four-lepton mass scale. With the 
current data there is no evidence for a statistically significant bias. 

The Crystal Ball width (sigma of the Gaussian core) returned by the fit is 1.1 ± 0.5 GeV, consis- 
tent with the simulation-based expectations of 1.0% for 4fi and 1.4% for lelji (which for a Higgs 
boson with mass 125 GeV corresponds to 1.3 GeV and 1.8 GeV) for a Higgs signal in this mass 
range. There are only two 4e events, so the mass resolution of 4e events cannot be constrained. 
With the current data, we can measure the average four-lepton mass resolution with about 45% 
statistical uncertainty. 

6 Summary 

The first observation of the Z boson decaying to four leptons (electrons and /or muons) in 
proton-proton collisions has been presented. The data set analyzed corresponds to an inte- 
grated luminosity of 5.02 fb -1 at s/s = 7TeV. A pronounced resonance peak, with a sta- 
tistical significance of 9.7 a, is observed in the distribution of the invariant mass of four lep- 
tons with mass and width consistent with expectations for Z boson decays. The event yields, 
branching fraction, and cross section reported here are defined by phase space restrictions on 
the leptons, namely, 80 < < 100 GeV, where is the invariant mass of the four lep- 
tons, and > 4 GeV for all pairs of leptons, where mu is the two-lepton invariant mass. 
We observe 28 events, in agreement with the expectation of 30.0 ± 2.6 events, comprised of 
29.4 ± 2.6 Z — > A£ events and 0.6 ± 0.2 events from backgrounds. The measured branch- 
ing fraction is B(Z U) = (4.2+{j| (stat. ) ± 0.2 (syst.)) x 10~ 6 , in agreement with the stan- 
dard model prediction of 4.45 x 10~ 6 . The measured cross section times branching fraction is 
cr(pp Z)B(Z -> U) = 112+|* (stat.)±5 (syst.)+2 (lumi.) fbr 1 , also consistent with the stan- 
dard model prediction of 120 fb. The four-lepton mass peak arising from Z — > A£ decays pro- 
vides a calibration channel for the Higgs boson search in the H — > ZZ — > A£ decay mode. 
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4e, Afi, and 2e2ji, are combined. The solid line represents a simultaneous fit to the background 
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